Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Pressure dependence of the resistivity around valence transition 
based on 1/iV-expansion study for extended Anderson lattice 

Yasutaka Nishida * and Tomokazu Ookubo t 

Department of Materials Engineering Science, 
Graduate School of Engineering Science, Osaka University, 
Toyonaka, Osaka 560-8531 

(Received February 5, 2008) 

We investigate the behavior of the resistivity around valence transition in the one-dimensional 
extended periodic Anderson model (PAM) with the Coulomb repulsion between f and conduc- 
tion electron Uf c . By using 1/iV-expansion in the leading order in (1/N)°, where N is the 
spin-orbital degeneracy of f-electrons, valence transition happens at critical Uf c with increasing 
the atomic f-level as seen in the previous work. 5 ' We calculate the resistivity whole tempera- 
ture range around valence transition based on the extended PAM including the crystal field and 
investigate how the physical properties such as Kondo temperature is modified by additional 
Ut c interaction. As a result, the double peak structure of the electrical resistivity fades away 
rapidly by the rapid increase of the Kondo temperature. 
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1. Introduction 

Heavy Fcrmion Systems have been investigated exten- 
sively both experimentally and theoretically as typical 
system of strongly correlated electrons. In 1979 Stcglich 
and collaborators discovered an intermetallic of the rare 
earth metal cerium CeCu2Si2, 1 - ) which is well known as 
first discovered superconductor with T c ~ 0.7K at am- 
bient pressure. In the phase diagram in P-T plane of 
CeCu2Si2 under the high pressure, two superconduct- 
ing transition temperature is confirmed. The isostructual 
compound, CeCu2Ge2, 2,3 ' also has the same phase dia- 
gram. 

These two compounds have been considered to have 
the similar physical properties if the origin of the pres- 
sure for CeCu 2 Si2 is shifted. According to the recent ex- 
perimental research, 4 ' the mechanism of the the super- 
conductivity at higher pressure region is seemed to be 
different from the ordinary antiferromagnetic spin fluc- 
tuation, because transition temperature is away from the 
point of magnetic instability at T = 0. Around the su- 
perconducting transition temperature at higher pressure 
P c . the residual resistivity has a peak, and the coeffi- 
cient of T 2 -tcrm of the resistivity shows rapid decrease. 
At high temperature, two temperatures corresponding to 
maxima of the magnetic resistivity coincides around P c . 
From these results, these coincidence around P c seemed 
to be related to a rapid change of valence of Ce, 5 ' 6 ^ so 
that the superconductivity mediated by the valence fluc- 
tuation was proposed. 

In this short paper, we study the pressure dependence 
of the resistivity based on the pioneering work 5 ) on the 
valence fluctuation by the extended periodic Anderson 
model with f-c coulomb repulsion Uf c . Then, we use 1/N- 
cxpansion 7-9 ' for extended Anderson lattice including 
crystal field, calculate both the low temperature region 
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(T < Tr-/10, where Tk is the Kondo temperature.) and 
the high temperature region (T > Tk) within (1/N)°. In 
the case of Uf c = 0, the double peak structure due to the 
crystal field is shown to merge into single peak gradually 
with increasing £/(< 0). On the other hand, in the case 
of U{ c 7^ 0, we find that the double peak structure fades 
away more rapidly due to the increase of Tk caused by 
rapid change of n/. 

The paper is organized as follows. In §2, we introduce 
the model Hamiltonian and l/A^-expansion. In § 3 and 4, 
we show the result the physical properties both low and 
high temperature region, and the pressure dependence 
of the resistivity. In section 5, we summarize conclusions 
and remarks. 

2. Model and Formal Preliminaries 

In this paper, we start with an extended periodic 
Anderson model, which is a periodic Anderson model 
(PAM) with f-electron in a manifold of J = 5/2, 10 > and 
Coulomb repulsion between f and conduction electrons, 
Uf c : Our model Hamiltonian is given by 

H = H c + H { + H hy + H {c , (2.1) 
ff c =]T£ k ^c ka , (2.2) 

H { = J2EirflfiT, (2.3) 
ir 

H hy = ^= Yl (Wr ik - R 'cL.Ar&l+h.c.), (2.4) 

V L a,i, k,r 

H lc = U lc J2<r<, (2-5) 

i,r 

where, c^ a is the creation operator for the conduction 
electron with wave vector k and spin a, and £kcr is 
the energy of the conduction electron with wave vec- 
tor k and spin a. Nt, denotes the number of the lat- 



2 J. Phys. Soc. Jpn. 



Full Paper 



Online- Journal Subcommittee 



tice site, and i stands for the site index. Furthermore 
we introduce the slave-boson which represents the / - 
state, and pseudo-fermion which represents /-"--state, fol- 
lowing the Coleman: b\ is the creation operator for the 
slave-boson at i-site, // r is the creation operator for the 
pseudo-fermion at i-site with Crystalline-Electric-Field 
(CEF) level The mixing potential Vkrv is expressed 
as Vkro- = Y^m OrAiVkMa, M is the z-component of an- 
gular momentum J of an f-electron: M — J Z (J = 5/2). 
Using the spherical harmonic function F ; m , VkMa- is given 
by 



Vi 



kA/o- 



Vq 



oVt{- 2 -i 



\l-2Mo) 



1 = 3 



(2.6) 



Otm is determined by the symmetry of the crystal field. 
In tetragonal symmetry we can put each f-level as follows, 



Ei\±l) — £ f + A2, Ei\± 2 ) 



T -f + Ai, Ei\±3) — Sf. 

(2.7) 



where, the CEF splitting is parameterized by A 12 - 

Here the on-site Coulomb repulsion U between f- 
electrons is assumed to be infinite, so that we must im- 
pose the following local constraints in order to guarantee 
an equivalence between the present model Hamiltonian, 



b]b< 



fir fir 



1. 



(2. 



The expectation value of an operator O under the local 
constraint eq. (2.8) is given as 11 -* 



(o)= hm (on^>A/<n^ 



(2.9) 



where (A)\ is calculated in the grand canonical ensemble 
for the Hamiltonian Hy. 



{A)x 
H x = 



= Tr[e-^l 



/Tr 



(2.10) 
(2.11) 



The single-particle Green function G^ a for the conduc- 
tion electron, Bi for the slave boson and Tir for the 
pseudo-fermion is given as follows: 

Gk CT (iw„) = (iw„ - e k(T - S k(T (iLj„)) _:L , (2.12) 

Bi(w n ) = (iv n - Xi - UiQvn))' 1 , (2.13) 

F ir (iuj n ) = (iu n — EiT — A,) -1 . (2.14) 

The self-energies Eko- and Ui is calculated by 1/N- 
cxpansion procedure. We perform an extended calcu- 
lation on the basis of the modified version of 1 /TV- 
expansion method. 9 ^ In the lowest order, the Dyson equa- 
tions for one-particle Green's function are given by the 
diagram shown in Fig. 1. Then the analytic form of the 
self-energy is given by 

Sk CT (iw n ) = hm -y^|Viiv| 2 2 , y^.Bi(ii/ n ) 

Ai— >oo L ^ — * *• — * 



x F° r (iuj n +iv n )/(Qi}\ 



(2.15) 




i f i,r \ i 



k ° 



B t : 
Fir ■ 

<Qi>% 



<Q<>x 



Fig. 1. Diagrammatic representation of Dyson equations for sin- 
gle Green function within accuracy of (l/N)° . Solid line denotes 
the conduction electron propagator, wavy line the slave-boson 
propagator and dashed line the pseudo-electron propagator, re- 
spectively. 




ct k„ 



x F° r (iu n + iu n ), 



(Qi) 



r 



JifiTjx + (n bi )\. 



where 



and 



(n fiT )\ = {fj T fir)x, 



{nu)\ 



(2.16) 
(2.17) 

(2.18) 
(2.19) 



By solving the coupled self-consistent equations eqs. 
(2.12)-(2.14) over the whole temperature range, we can 
discuss the temperature dependence of the physical 
quantities. In the next section, we discuss the physical 
properties at T = as first step. 

3. Physical properties at low temperature limit 

At lower temperatures than the coherence temperature 
T , the coupled self-consistent equations become more 
simplified, with the help of the following relation, 

--ImBi(u + i0+) ~ ad(uj - X. t - e f + E ) + C(uj), 

(3.1) 

where, Eq and a are the binding energy and residue of the 
slave-boson, respectively, and C{lo) is continuum part, 
which is non-zero oj > Aj + £/. Eq corresponds to the 
Kondo temperature Tk, and the coherence temperature 
To is estimated with Eq/10. Eq and a are determined by 
the relation as 



E =e f - ReILi(Xi +e f -Ea) 



= 1 



(RelliM) 



L0 = \i+ef-E 



(3.2) 
(3.3) 



Here we treat £/ fc term with the Hartree-Fock approx- 
imation in O((l/A0 ). 12 ) By using eq. (3.1), the system 
is in the coherent regime with the renormalized bands 
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given by 



Gko-(i^n) = E ' 



(3.4) 



*kcr 



with 

£k CT = £k CT + Ut c rif, (3.5) 
= £W + ReE k(T (a k(T ), (3.6) 

(«L - #o)(aL - - Ai)("L - - A 2 ) 



(4 CT -4^)(4 CT -4f)(4 CT -4t 3 ) 



(3.7) 



where, a k(7 = a k + 4 . 

Hence the self-consistent equations cqs. (2.12)-(2.14) 
are written by 



^LlW/(6L) 



(3.8) 



^ = 1 +^-EEEE 



j=l r a- k„ 



(£ 0+J E; jr - e/ -a k J 



(3.9) 



= ^ E E E /(«L)^L + (i - «)■ 



j'=i <y k„ 



and 



-E^r = Eir + U{ c n c , 



(3.10) 
(3.H) 



where n, n c , and rt/ is the total electron number, the 
number of the conduction electrons, and f-electrons per 
site, respectively. Finally, the binding energy Eq, the 
residue a and other physical properties at T = are 
obtained from a series of equations eqs. (3.8)-(3.10). 
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Fig. 2. (a) rif vs sj in the case of n = 1.4, Uf c = and V 2 = 
0.98 x 10" 2 D 2 for a series of CEF splittings, (b) rif vs Ef in the 
case of n = 1.4, U{ c = and V 2 = 2.8 X 10 _2 D 2 for a series of 
CEF splittings. 
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3.1 f-electron number 

In this section, we show the f-electron number n f per 
site as a function of the atomic f- level e / for series of CEF 
splitting scheme with Uf c . At first, let us investigate the 
trivial case, V 2 = J2* l^kivl 2 = 0.98 x lO" 2 !? 2 , n = 
n c + nf = 1.4, and Ai = A 2 = 0.6-D. Then the grand 
state is doublet (Acs = 2) state, so that, in Fig. 2, we 
find that if Ef < —0.6, rif ~ 1. However in case of V 2 = 
2.8x 10~ 2 , rif decreases monotonously by increasing the 
width of the density of state of localized f-electrons due 
to hybridization effect. On the other hand, in the case 
of Ai = A 2 = (Acs =6), nf decreases monotonously. 
If TVgs is small, rif shows a little sudden decrease with 
increasing e / . 

Next, we show the Uf c ^ result in Fig. 3 and 4. In 
Fig. 3, we set Ai = A 2 = for simplicity. Then, rif 
shows the rapid decrease gradually with increasing C/ fc . 
This behavior was already obtained by the past pioneer- 
ing work. 5 ) In Fig. 4, we set Ai = 0, A 2 = 0.3D and 
Ai = A 2 = 0.6Z? in order to consider the CEF effect. 
From these result we find rif vs £/ is not much modified 
qualitatively by CEF effect. 



Fig. 3. rif vs Sf in the case of n = 1.4, V 2 = 2.0 X 10 2 D 2 and 
Ai = A2 = for a series of C/f c . 

3.2 Kondo temperature 

The Kondo temperature corresponds to the binding 
energy of the slave-boson Eq in lattice case. We present 
the solution Eq for the series of the self-consistent equa- 
tions (3.8)-(3.10) within the 0((1/A)°). Incidentally Eq 
can be written as follows in the case of impurity without 
CEF effect and U lc , 

* = ^(^V*)' (3 - 12) 

where we impose the condition that D S> |e/| 3> E . 
In eq. (3.12), 6 means the degeneracy of f-electron in a 
manifold J = 5/2, and po represents for the density of 
state of the conduction electron per spin. 

In Fig. 5, we show the result of Eq vs £/ for the vari- 
ous Uic and CEF splitting. At first we again consider the 
trivial case Ai = A 2 = and U{ c = 0. Then, Eq increase 
gradually with increasing £/(< 0). This tendency is ex- 
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degeneracy of the localized f-electrons. 




Fig. 4. (a) rij vs Ef in the case of n = 1.4, V 2 = 2.0 X 10 _2 D 2 , 
Ai = 0, and A2 = 0.3D for a series of U{ c . (b) rif vs Ef in the 
case of n = 1.4, V 2 = 2.0 x 10~ 2 D 2 , Ai = 0.6£>, and A 2 = 0.6D 
for a series of Uf c . 



Fig. 6. rif vs Ef and i?o vs in the case of n = 1.4, V = 
2.0 X 10" 2 D 2 . The result of Ai = A 2 = 0.6D is shown by solid 
line, and that of Ai = A2 = is shown by enclosed circle. 



pcctcd from the impurity result eq. (3.12). Next, in case 
of Uf c 7^ 0, Eq changes drastically as large as n; with in- 
creasing e f . This is explained by the enhancement of the 
renormalization factor q accompanying with the rapid va- 
lence change. The renormalization factor q derived from 
mean-field 13 ^ and Variational Monte Cairo 14 - 1 is obtained 
as follows 



g _x = 1 ~ n//2 



(3.13) 

1-n/ 

The rapid change rif : 1 — *■ causes also the rapid en- 
hancement of q, so that the Kondo temperature is en- 
hanced around valence transition. 
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Fig. 5. E vs e f in the case of n = 1.4, V 2 = 2.0 X 10 _2 D 2 , 
Ai = A2 = for a series of Cf c . 



In Fig. 6, we illustrate the result of the Ai = A2 = 0.0 
and Ai = A2 = 0.6D under the condition Uf c = 0. In 
the previous subsection, we find the CEF effect does not 
affect the £/ dependence of rif, on the other hand, Kondo 
temperature Eo is greatly influenced of CEF effect. The 
value of the Kondo temperature depends much by the 



4. Electrical Resistivity 

In this section, the effect of the pressure on the tem- 
perature dependence of the resistivity is discussed. The 
pressure effect is parameterized as the atomic f-level £/ 
measured from the band center of conduction electrons. 
The conductivity is obtained by means of the Kubo for- 
mula. By omitting some constant factors, we define the 
reduced conductivity, 



1 \ ^ 2 

er,k 



de - 



dm 

de 



ImGkcr (e + i0-| 



(4.1) 



where we use the renormalized Green function derived 
from 1/iV-expansion procedure. In order to calculate the 
resistivity 1/a over the whole temperature range, we 
must solve the Dyson equation eqs. (2-12)-(2.14) self- 
consistently (self-consistent 1/N). 

We illustrate the result of the calculation for the series 
of ef in Fig. 7 The parameters are adopted as £>=l(band 
width), V 2 = J2a l^kral 2 = 1.28xlO- 2 L> 2 , n = n c +n f = 
1.4, and set the CEF parameter as At = 0.015D and 
A 2 = 0.038D putting the case of CeCu 2 Si 2 . 15) 

In both case of Uf c = and Uf c ^ 0, the double peak 
structure fades away with increasing £/ (< 0) . The double 
peak structure of the temperature dependence is com- 
monly shown to merge into a single peak with increasing 
pressure. These behaviors are consistent with the ten- 
dency observed in Ce-based heavy fermions. However, at 
the large £f (high pressure region), the position of the 
single peak of 1/a differs much in both cases. In case 
of U{ c = 0, the peak position does not change much at 
large £f. On the other hand, in case of U& 7^ 0, the 
peak position changes much with increasing Ef. The Fig. 
9 shows this situation clearly, i.e., rapid increase of Eq 
(Kondo temperature) happens around the valence tran- 
sition. Under the condition without Ufc, the Kondo tern- 
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Fig. 9. nf vs £f and Eo vs £/ in the case of n = 1.4, V = 
1.28 X lCr 2 D 2 , Ai = 0.015D and A 2 = 0.038D for the series of 
U tc . 



Fig. 7. (a) The temperature dependence of electrical resistivity 
for a series of the atomic f-level in the case of U{ c = and (b) 
that in the case of Uf c = 0AD. The CEF splittings are set as 
Ai = 0.015D and A 2 = 0.038D in both cases. 
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Fig. 8. The temperature dependence of electrical resistivity for a 
series of the atomic f-level in the case of U{ c = 0.5D. 



perature increases slowly with increasing Sf compared 
with the case of Uf c = 0AD. The effect of the additional 
term Uf c seems to promote the merging of the double 
peak structure of the resistivity. 

We also illustrate the result in case of U[ c = 0.5D in 
Fig. 8. We find that the electrical resistivity becomes 
more sensitive to the change of £/ than that of Uf c = 0. 

5. Conclusion 

We have studied the behavior of the pressure depen- 
dence of the electrical resistivity around valence transi- 
tion based on 1 /TV-expansion for extended periodic An- 
derson lattice including CEF effect. This is the extended 
calculation of the pioneering work. 5 ) At lower temper- 
ature than Kondo temperature, the result of nj vs £/ 
does not change qualitatively, i.e., nf shows the rapid 
decrease with increasing Uf c where CEF effect does not 



affect the tendency of the change of n/ as a function of Ef . 
remarkably. However the Kondo temperature Tk shows 
different behavior by CEF effect, because Tk depends on 
the degeneracy of the localized f-electrons much. More- 
over, we have calculated the electrical resistivity over the 
whole temperature region by solving the Dyson equation 
self-consistently within (1/N)° and investigated the ef- 
fect of the additional interaction U{ c to the temperature 
dependence of the electrical resistivity. 

As a result, we found that the double peak structure 
of the resistivity due to the crystal field fades away more 
rapidly caused by the rapid change of n/. Especially, in 
the region of Sf around the valence transition, the resis- 
tivity show the drastic change by a minute change of £/ . 
The effect of the additional term Uf c promote the merg- 
ing of the double peak structure of the resistivity This 
tendency may be consistent with the merge of the resis- 
tivity right above P c in the phase diagram in P-T plane 
of CeCu2Si2 or CeCu2Ge 2 . 

Finally we point out some remaining problem and per- 
spectives. In our calculation, the term of U[ c is treated 
within Hartrce-Fock approximation, so that the effect of 
Ui c may be overestimated by the present calculation. In 
order to discuss the problem more qualitatively, we need 
more proper approximation. Experimentally, the rapid 
decrease of Kadowaki- Woods ratio is observed around 
P c . This suggests that the system is changed from Kondo 
regime to the valence fluctuation regime rapidly. We can 
discuss this problem by studying up to (1/7Y) 1 , where 
intersite correlation can be discussed. By investigating 
between the change of the degeneracy of the f-electron 
and pressure effect or Uf c , we may also discuss the mech- 
anisms of enhancement of Tk in detail. 
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